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ABSTRACT 

The scope of this paper is to compare the catalog extraction performances obtained using 
the new combination of SExtractor with PSFEx, against the more traditional and diffuse 
application of DAOPHOT with ALLSTAR; therefore, the paper may provide a guide for the 
selection of the most suitable catalog extraction software. Both software packages were tested 
on two kinds of simulated images having, respectively, a uniform spatial distribution of sources 
and an overdensity in the center. In both cases, SExtractor is able to generate a deeper catalog 
than DAOPHOT. Moreover, the use of neural networks for object classification plus the novel 
SPREAD_MODEL parameter push down to the limiting magnitude the possibility of star/galaxy sep- 
aration. DAOPHOT and ALLSTAR provide an optimal solution for point-source photometry in 
stellar fields and very accurate and reliable PSF photometry, with robust star-galaxy separation. 
However, they are not useful for galaxy characterization, and do not generate catalogs that are 
very complete for faint sources. On the other hand, SExtractor, along with the new capability 
to derive PSF photometry, turns to be competitive and returns accurate photometry also for 
galaxies. We can assess that the new version of SExtractor, used in conjunction with PSFEx, 
represents a very powerful software package for source extraction with performances comparable 
to those of DAOPHOT. Finally, by comparing the results obtained in the case of a uniform 
and of an overdense spatial distribution of stars, we notice, for both software packages, a de- 
cline for the latter case in the quality of the results produced in terms of magnitudes and centroids. 
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1. Introduction 

Over the past two decades, advances in tech- 
nology are progressively moving us beyond the 
traditional observational paradigm, in which most 
astronomical studies were made by individual ob- 
servations of small samples of objects. Modern 
digital sky surveys are already producing increas- 
ing amounts of data flows that cannot be effec- 
tively handled with traditional methods. This 
has profoundly changed the needs of scientists in 
terms of software and data-analysis methods. For 
instance, the sheer size of the raw data makes it 
almost impossible to re-process the raw images 
and, therefore, large catalogs are becoming the 
primary source of information. 
When extracting a catalog, the main aspects to 
take into account are: to detect as many sources 
as possible; to minimize the contribution of spuri- 
ous objects; to correctly separate sources in their 
classes (e.g. star/galaxy classification); to produce 
accurate measurements of photometric quantities; 
and, finally, to obtain accurate estimates of the 
positions of the centroids of the sources. 
Among the main source extraction software pack- 
ages used by the astronomical community, there 
are SExtractor (Bcrtin fc Arnoutsl Il99q) and 
DAOPHOT II (|Stetsonlll987t ). and the latter is 
often used in co mbination wit h its companion 
tool ALLSTAR (jStetsonl Il994l ). SExtractor is 
commonly used in extragalactic astronomy and 
has been designed to extract a list of measured 
properties from images, for both stars and galax- 
ies. DAOPHOT and ALLSTAR were designed to 
perform mainly stellar photometry. So far, only 
DAOPHOT II, used together with ALLSTAR, has 
been able to produce more accurate photometry 
for stellar objects using the Point Spread Function 
(PSF) fitting technique, while the PSF fitting pho- 
tometry in SExtractor has, instead, become pos- 
sible only in the recent years. The first attempt 
was in the late 90s, when the PSFEx (PSF Ex- 
tractor) software package became available within 
the TERAPIX "consortium". This tool extracts 
precise models of the PSF from images processed 
by SExtractor. Only after 2010, th rough the pub- 
lic release of PSFEx (|BertirJl201 lH 1 !. and with the 
recent evolution of computing power, has PSF 
fitting photometry become fully available in SEx- 



tractor. 

The scope of this paper is to compare the re- 
sults obtained using the combination of SEx- 
tractor with PSFEx, and DAOPHOT with ALL- 
STAR, by focusing, in particular, on the com- 
pleteness and reliability of the extracted cata- 
logs, on the accuracy of photometry, and on the 
determination of centroids, both with aperture 
and PSF-fitting photometry. A previous com- 
parison a mong extraction softw are tools was per- 
formed by iBecker et all (|2007h . They, in pursuit 
of LSST science requirements, performed a com- 
parison among DAOPHOT, two versions of SEx- 
tractor (respectively 2.3.2 and 2.4.4), and DoPhot 
(|Mateo fc Schechteij[l989l) . However, differently 
from the present work where simulations are used, 
they evaluated as "true" values the measurements 
obtained with the S PSS imaging pipeline photo 
([Lupton et al.1 120011 ). Furthermore, we wish to 
stress that their results were biased by the fact 
that in 2007 the PSF fitting feature had not yet 
been implemented in SExtractor. 
The present work performs, for the first time, 
a comparison between DAOPHOT and SExtrac- 
tor PSF photometry, providing a guide for the 
selection of the most suitable catalog extraction 
software packages. 

The simulations used for the comparison are de- 
scribed in Sect. [2] In Sect. El the main input pa- 
rameters of the software packages are overviewed 
and the adopted values are specified in Sect.|U In 
Sect. the obtained results are shown. In order 
to better evaluate the performances of both soft- 
ware packages on crowded fields, we describe: in 
Sect.[6j a test performed on an image showing an 
overdensity in the center. Finally, the results are 
summarized in Sect. [JJ together with our conclu- 
sions. 



2. Image Simulations 

Image simulations are suitable in testing perfor- 
mances of various software packages. Simulations, 
in fact, allow us to know exactly the percentage 
and the type of input sources and their photomet- 
ric properties. 

In the present work, simulations have been ob- 
tained by using two software packages: Stufg and 



1 Available at http: //www. astromatic .net/sof tware/psf ex 



2 Available at http: //www. astromatic .net/sof tware/stuff 
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SkyMakeH developed by E. Bertin. With these 
tools, it is possible to reproduce the real outcome 
of a CCD observation, once the characteristics 
of the telescope and the camera are known. In 
practice, Stuff can be used to produce a realistic 
simulated galaxy catalog, while SkyMaker uses 
this catalog to produce an optical image under re- 
alistic observing conditions, allowing also to add 
a stellar field. 

The galaxy catalog simulated by Stuff can be pro- 
duced so to be consistent with the assumed cosmo- 
logical model and with the statistical distributions 
of stars and galaxies in terms of redshift, luminos- 
ity and color. In a binned redshift space, Stuff 
produces galaxies of different Hubble types E, SO, 
Sab, Sbc, Scd and Sdm/Irr. The number of galax- 
ies in each bin is determined from a Poisson dis- 
tribution, by assumi ng a non-evolvi ng Schechter 
luminosity function ([Schechter 1976). Cosmolog- 
ical parameters, luminosity function, as well as 
instrumental parameters are specified by the user 
in the input configuration file. In particular, the 
size of the image, the pixel scale, the detector gain, 
and the observed passbands must be provided as 
input information. Filters can be selected among 
the many available in the wavelength range [0.29, 
87.74831] /an. Finally the magnitude range of 
simulated galaxies has to be fixed. 
The image is then created by SkyMaker, by ren- 
dering sources of the input catalog in the frame at 
the specified pixel coordinates, to which is added 
a uniform sky background, Poissonian noise, and 
Gaussian read-out noise. Stellar sources are mod- 
eled using a PSF internally generated, while the 
various types of galaxies are modeled as differently 
weighted sums of a bulge profile and an exponen- 
tial disk. The PSF profile takes into account both 
atmospheric seeing and optical aberrations. There 
are many parameters to be set in the SkyMaker 
configuration file. The most important are related 
to: the pupil features, e.g., the size of the mirrors 
and the aberration coefficients; the detector char- 
acteristics, e.g., gain, saturation level and image 
size; and the observing conditions, e.g, full width 
at half maximum of the seeing (FWHM) and ex- 
posure time. 

For our work, we simulated images as they would 



be observed by the VSlfl (VLT Survey Telescope) 
and the OmegaCAM camer£0. The field of view 
(FoV) of OmegaCAM@VST is 1 square degree, 
with a pixel scale of 0.213 arcsec/pixel. In order 
to reduce the computational time, we limited our 
simulations to a FoV of 1/4 of VST. The aberra- 
tion coefficients, the tracking errors, and the po- 
sitions of the spiders were set properly, according 
to the VST technical specifications. The allowed 
magnitude range was set to 14-26 mag and the 
exposure time was fixed at 1500 s. Finally, ac- 
cording to the ESO statistics at Cerro Paranal, 
the FWHM of the seeing was set to 0.7 arcsec. In 
this paper, we report results relative to B-band 
simulated images. We obtain a catalog of N=4120 
sources down to the input magnitude limit. This 
corresponds to an image with an average surface 
number density of ~ 20 sources/arcmin 2 . 
In order to perform a more complete comparison 
of both software packages, we also simulated a 
crowded image with an overdensity in the center 
(see Sect. [6]). To this purpose, the parameters of 
Stuff and SkyMaker are equal to those used for 
images with uniform source distribution, except 
for the image size and the exposure time. For 
these second group of tests, we reduced the expo- 
sure time to 300 s, with an input magnitude limit 
of 14-25, while the image size is set to 2048x2048 
pixels (~ 7x7 arcmin 2 ). We wish to stress that, 
in the simulations, we did not include any arti- 
fact such as bad pixels, ghosts or bad columns 
or other effects that, while being crucial in other 
applications, are of no interest here, since we are 
comparing the performances of two different soft- 
ware packages. 



3. Source extraction software 

In the following section, we briefly discuss 
how DAOPHOT works in combination with ALL- 
STAR, and how SExtractor works in combination 
with PSFEx, and give a brief overview of the main 
parameters that are needed to be set in order to 
optimally run the selected software packages. 



4 http: //www. eso . org/public/teles-instr/surveytelescopes/vst /surveys .html 
5 http: //www. astro- wise . org/~omegacam/index . shtml 



3 Available at http: //www. astromatic .net/sof tware/skymaker 
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3.1. DAOPHOT II 

DAOPHOT II is composed of a set of routines 
mainly designed to perform stellar photometry 
and astrometry in crowded fields. 
It requires several input parameters, listed in the 
file daophot.opt, including detector gain, readout 
noise (GAIN, READ NOISE), saturation level (HIGH 
GOOD DATUM), approximate size of unresolved 
stellar sources in the frame (FITTING RADIUS), 
PSF radius (PSF RADIUS), PSF model (ANALYTIC 
MODEL PSF), and a parameter designed to allow 
the user to visually inspect the output of each 
routine (WATCH PROGRESS). 

The first step performed by DAOPHOT II is to es- 
timate the sky background and to find the sources 
above a fixed input threshold through the FIND 
routine. 

This threshold represents the level (in ADU), 
above the sky background, required for a source to 
be detected. In order to ignore smooth large-scale 
variations in the background level of the frame, 
the image is convolved with a lowered truncated 
Gaussian function, whose FWHM is equal to the 
input value set by the FWHM parameter. After the 
convolution, the program searches for the local 
maxima sky enhancement. 

Once the sources are detected, DAOPHOT II per- 
forms aperture photometry via the PHOTO rou- 
tine. Aperture photometry usually requires the 
definition of at least two apertures. The first one 
is usually circular, centered on the source and with 
a radius of a few times its FWHM. The second one 
is, instead, ring-shaped; usually this is concentric 
to the first one and it has an inner radius equal to 
the radius of the first aperture. The ring-shaped 
aperture is used to estimate the sky contribution 
and it usually covers a number of pixels equal to 
or at least comparable with that of the aperture. 
Then, the flux of the source is obtained by sub- 
tracting the sky flux from the aperture flux. The 
aperture size must be chosen thoroughly. In fact, 
if the radius of the inner aperture is too small, 
there will be a flux loss; while, if it is too large, 
too much sky is included and the measurements 
will become too noisy. The radii of the aperture 
and sky annulus for DAOPHOT can be specified 
in an input file: photo, opt. The inner and outer 
radii of the sky annulus, which is centered on the 
position of each star, must also be specified. 
Beside the source magnitude, the PHOTO routine 



produces the coordinates of the source centroids, 
corresponding to the barycenters of the intensity 
profile around the source. 

Aperture photometry performs rather well un- 
der the hypothesis of bright and isolated stars. 
However, the stars in crowded fields are faint and 
tend to overlap. In these cases, the PSF fitting 
photometry can produce better results. This mea- 
surement requires a PSF model to be derived from 
the stars in the image of interest. The normalized 
PSF model is then fitted to each star in the image, 
in order to obtain the intensity and magnitude. 
DAOPHOT II can build a PSF model from a 
sample of stars obtained with the PHOTO rou- 
tine in an iterative procedure, intended to subtract 
neighboring stars that might contaminate the pro- 
file. Among them, DAOPHOT will exclude stars 
within one radius from the edges of the image and 
the stars too close to saturated stars. The ana- 
lytical formula of the PSF is chosen by the user 
among several available models: a Gaussian func- 
tion, two implementations of a Moffat function, 
a Lorentz function , and two im plementations of 
a Penny function (jPennvl [l995h . The PSF rou- 
tine produces a PSF model and a list of the PSF 
stars and their neighbors. The modeled PSF stars 
can be visually inspected by setting properly the 
aforementioned WATCH PROGRESS parameter. 
Although DAOPHOT is designed for stellar 
photometry, extended sources are likely always 
present in real images and, therefore, a reliable 
method is required to separate galaxies from stars. 
The sharpness parameter (SHARP) can be used as 
star/galaxy classifier. SHARP describes how much 
broader is the actual profile of the object com- 
pared to the profile of the PSF. The sharpness 
is, therefore, dependent on the model of the PSF, 
and can be easily interpreted by plotting it as 
a function of apparent magnitude. Objects with 
SHARP significantly greater than zero are probably 
galaxies. 

Throughout the paper, we indicate simply with 
DAOPHOT the stand-alone DAOPHOT II ver- 
sion 

3.2. ALLSTAR 

After having derived PSF models with DAOPHOT, 
ALLSTAR simultaneously fits multiple overlap- 



3 Available at http://starlink.jach.hawaii.edu/starlink 
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ping point-spread functions to all the detected 
sources in the image. At each iteration, ALLSTAR 
subtracts all the stars from a working copy of the 
input image, according to the current best guesses 
about their positions and magnitudes. Then, it 
computes the increments to the positions and 
magnitudes by examining the subtraction resid- 
uals around each position. Finally, it checks each 
star to sec whether it has converged or it has be- 
come insignificant. When a star has converged, its 
coordinates and magnitude are written in the out- 
put file, and the star is permanently subtracted 
from the working copy of the image; when a star 
has disappeared, it is simply discarded. 
The input parameters for ALLSTAR, listed in all- 
star, opt, are similar to those in daophot.opt and 
photo, opt. 

Moreover, it is possible to optimize the determi- 
nation of centroids, by applying a PSF correction 
and setting the option REDETERMINE CENTROIDS. 
An improvement of the star/galaxy classification 
is possible by using the aforementioned sharpness 
measure obtained by ALLSTAR (SHARP). This 
value may be used also in conjunction with an- 
other ALLSTAR output parameter, for instance 
X, which is the observed pixel-to-pixel scatter 
from the model image profile, divided by the ex- 
pected pixel-to-pixel scatter from the image pro- 
file. Throughout the paper, we indicate simply 
with ALLSTAR the PSF-fitting software package 
that comes together with DAOPHOT II version 
v. 1.3-2. 

3.3. SExtractor 

SExtractor is a software package mainly de- 
signed to produce photometric catalogs for a large 
number of sources, both point-like and extended. 
Sources are detected in four steps: i) sky back- 
ground modeling and subtracting, ii) image fil- 
tering, Hi) thresholding and image segmentation, 
iv) merging and/or splitting of detections. The 
final catalog is, indeed, extracted according to the 
input configuration file, in which parameters are 
set by the user. 

The first step of the background estimation can 
be skipped if the user provides manually an input 
estimation of sky background. For the automatic 
background estimation, the most critical input pa- 
rameters to be set are BACK_SIZE, the size of each 
mesh of the grid used for the determination of 



the background map, and BACK FILTERSIZE, the 
smoothing factor of the background map. 
Once the sky background is subtracted, the im- 
age must be filtered. This implies convolving the 
signal with a mask, shaped according to the char- 
acteristics that the user wants to enhance in the 
image data. In fact, there are different filters 
available in SExtractor. The more suitable are 
"top-hat" functions, optimized to detect extended, 
low-surface brightness objects, Gaussian functions 
usually used for faint object detection, and "Mex- 
hat" filters, which work with a high value of de- 
tection threshold, suitable for bright detections in 
very crowded star fields. 

The detection process is mostly controlled by the 
thresholding parameters (DETECT_THRESHOLD and 
ANALYS I S_THRESHOLD) . The choice of the thresh- 
old must be carefully considered. A too high 
threshold determines the loss of a high number 
of sources in the extracted catalog, while a too 
low value leads to the detection of spurious ob- 
jects. Hence it is necessary to reach a compromise 
by setting these parameters according to the im- 
age characteristics, the background RMS, and also 
to the final scientific goal of the analysis. 
Two or more very close objects can be detected as 
a unique connected region of pixels above thresh- 
old and, in order to correct for this effect, SEx- 
tractor adopts a deblending method based on a 
multi-thresholding process. Each extracted set 
of connected pixels is re-thresholded at N lev- 
els, linearly or exponentially spaced between the 
initial extraction threshold and the peak value. 
Also, here a compromise is needed to be found, 
since a too low value for the deblending parameter 
leads to a lack of separation between close sources, 
while a too high value leads to split extended faint 
sources in more components. Alternatively, it is 
possible to extract the catalog with different de- 
blending parameters and to merge detections for 
extended sources or close pairs. 
Once sources have been detected and deblended, 
the software tool starts the measurement phase. 
SExtractor can produce measurements of posi- 
tion, geometry, and of several types of photomet- 
ric parameters, including different types of magni- 
tudes. Among photometric quantities, there are: 
the aperture magnitude (MAG_APER), having the 
same meaning as explained in Sect. 13.11 the Kron 
magnitude (MAG_AUTD), which is the magnitude 
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estim ated through an adaptive aperture ( Kronl 
1980), and the isophotal magnitude (MAG_IS0), 
computed by considering the threshold value as 
the lowest isophote. 

Among the available position parameters, it is 
important to mention the barycenter coordinates 
(X_IMAGE, Y_IMAGE), computed as the first order 
moments of the intensity profile of the image, 
and windowed positional parameters (XWIN_IMAGE, 
YWIN_IMAGE), computed in the same way as the 
barycenter coordinates, except that the pixel val- 
ues are integrated within a circular Gaussian win- 
dow as opposed to the object's isophotal footprint. 
To separate extended and point-like sources, it is 
possible to use the stellarity index (CLASSJ3TAR), 
which results from a supervised neural network 
that is trained to perform a star/galaxy classifi- 
cation. CLASS_STAR can assume values between 
and 1. In theory, SExtractor considers objects 
with CLASS_STAR equal to zero to be galaxies, and 
those with value 1 as a star. In practice, stars are 
classified by selecting a CLASS_STAR value above 
0.9. Two other parameters, often used to dis- 
criminate between star and galaxies, are the half- 
light radius (FLUX_RADIUS), and the peak surface 
brightness above background (/i ma x)- When plot- 
ted against the Kron magnitude, these two pa- 
rameters identify a so-called stellar locus. 
Throughout the paper, we indicate simply with 
SExtractor the software version 2.14.7 (trunk. r284). 

3.4. PSFEx 

The last version of SExtractor can work in com- 
bination with PSFEx, a package able to build a 
model of the image PSF. The latter is expressed 
as a sum of NxN pixel components, each one 
weighted by the appr opriate factor in t he poly- 
nomial expansion (see iMohr et al.1 120121 ). Then, 
SExtractor takes the PSFEx models as input and 
uses them to carry out the PSF-corrected model 
fitting photometry for all sources in the image. 
PSFEx requires, as input, a catalog produced by 
SExtractor to build a model of the image PSF, 
which can be read back in a second run by SEx- 
tractor itself. In order to allow PSFEx to work, 
the first catalog produced by SExtractor must 
contain at least a given number of parameters, as 
explained in the PSFEx manuaj^. In particular, 



the catalog must contain the parameter VIGNET, 
a small stamp centered on each extracted source, 
used to model the PSF. The size of VIGNET must 
be chosen accordingly to the size of the photomet- 
ric apertures defined by PHOT_APERTURES . 
PSFEx models the PSF as a linear combination of 
basis vectors. These may be the pixel basis, the 
Gauss-Laguerre or Karhunen-Loeve bases derived 
from a set of actual point-source images, or any 
other user-provided basis. The size of the PSF 
and the number and type of the basis must be 
specified in the configuration file. 
By using SExtractor combined with PSFEx, it 
is possible to obtain various estimates of the 
magnitude, in addition to those described in 
the previous section: the magnitude result- 
ing from the PSF fitting (MAG_PSF), the point 
source total magnitude obtained from fitting 
(MAG_P0INTS0URCE), the spheroidal component of 
the fitting (MAG_SPHER0ID), the disk component of 
the fitting (MAGJDISK), and the sum of the spheroid 
and disk components (MAG_M0DEL). Moreover, it is 
also possible to measure morphological parameters 
of the galaxies, such as spheroid effective radius, 
disk aspect ratio, and disk-scale length. 
The model of the PSF may be employed to ex- 
tract a more accurate star / galaxy classification us- 
ing the new SExtractor classifier, SPREAD_M0DEL, 
which is a normalized, simplified linear discrimi- 
nant between the best fitting local PSF model and 
a more extended model made by the same PSF 
convolved with a circular exponential disk model 
with scalelength = FWHM/16, where FWHM is 
the full- width at h alf maximum of the PSF model 
(jDesai et al.ll2012j ). 



A more detailed description of PSFEx a nd the 



new S Extr actor capabilities can be found in lBertin 
(|201ll) and lArmstrong et all (|2010f ). 
Throughout the paper, we indicate simply with 
PSFEx the software version 3.9.1. 

4. Catalog extraction 

In this section we provide a general discussion 
on how to set the input parameters in order to 
extract the catalogs with the software tools pre- 
sented above. 



'https : //www. astromatic .net/pubsvn/sof tware/psf ex/trunk/doc/psf ex.pdf 
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4.1. DAOPHOT and ALLSTAR 

Besides instrumental parameters, such as gain, 
saturation level, and readout noise, which are set 
according to the values used for the simulations 
(Sect. HJ, in DAOPHOT the detection and photo- 
metric options must be configured by means of in- 
put setup files. In the present analysis, the thresh- 
old value was chosen to detect as many possible 
sources, while avoiding as much as possible spuri- 
ous detections. In fact, as the threshold decreases, 
the number of detected sources increases up to a 
certain value, for which the relation changes in 
steepness. Thus it is possible to choose a rea- 
sonable value for the threshold, by plotting the 
number of extracted sources for different thresh- 
old values and by choosing the threshold near the 
"elbow" of the function. Moreover, in order to 
avoid spurious detections, due to Poissonian noise, 
the extracted catalog was visually inspected. In 
fact, by using the input parameters as reported in 
Tab[IJ only the 5% of the objects extracted in the 
catalog is spurious. By decreasing the threshold, 
although the number of extracted input sources 
increases, the percentage of spurious detections 
rapidly grows up to 30% at 4er threshold, even 
reaching more then 88% at 3a. 
The FITTING RADIUS was set equal to the FWHM 
of the image (see Tab. [I]). 

To obtain the best aperture radius, we have de- 
rived the growth curve for input stellar sources. 
Then, we fixed the aperture radius to 12.5 pixel 
(see al in Tab. [T]), producing the better cover- 
age of the sources input magnitude. Thus, values 
of INNER RADIUS and OUTER RADIUS were chosen 
accordingly, smaller and greater than the aperture 
radius, respectively. 

The PSF analytical model was chosen with the 
higher level of complexity, that is, the implemen- 
tations of the Penny function with five free param- 
eters (option 6). We chose to visually inspect the 
image of the PSF produced by DAOPHOT for all 
the PSF stars. ALLSTAR parameters were set ac- 
cordingly to those established for DAOPHOT and, 
furthermore, we required the redetermination of 
the centroids. 

The main parameters set for DAOPHOT and 
ALLSTAR are reported in Tab. HJ 



Parameter 


Values 


FITTING RADIUS 


3.38 


THRESHOLD (in sigmas) 


5 


ANALYTIC MODEL PSF 


6 


PSF RADIUS 


7.5 


al 


12.5 


INNER RADIUS 


20 


OUTER RADIUS 


35 


REDETERMINE CENTROIDS 


1.00 



Table 1. Main input parameters set in the 
DAOPHOT and ALLSTAR configuration files. 
FITTING RADIUS, PSF RADIUS, al, INNER RADIUS 
and OUTER RADIUS are expressed in pixels. 

4.2. SExtractor and PSFEx 

As for DAOPHOT, SExtractor instrumental 
parameters are set accordingly to those defined as 
input in the simulations (see Sect. [2]). 
Concerning the sky background modeling and sub- 
traction, we decided to automatically estimate the 
background within the software package adopting 
the global background map. Given the average 
size of the objects, in pixels, in our images, we 
chose to leave BACK_SIZE to the default value 64. 
The choice of the filter was more complex. We 
performed several tests with various filters ob- 
taining the best results by Gaussian and top-hat 
masks. However, the choice between the various 
filters, although affecting the number of detected 
sources, does not alter their measurements. 
For the thresholding parameters, we followed the 
same procedure approached with DAOPHOT, as 
described in Sect. 14. 1[ that is, by choosing a value 
near to the change in gradient of the relation be- 
tween the number of extracted sources and the 
threshold value for detections. Moreover, the 
catalog was visually inspected to avoid residual 
spurious detections and to verify the deblending 
parameters. 

We fixed the size of the aperture for photome- 
try, according to the one set in DAOPHOT, to 25 
pixels of diameter (PHOT_APERTURES). For PSFEx 
parameters we used a set of 20 pixel basis and a 
size for the PSF image of 25 pixels according with 
the aperture size. We adopted a 25x25 pixel ker- 
nel following PSF variations within the image up 
to 2 nd order. The main values set for SExtractor 
and PSFEx are reported in Tab. [2] 
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Parameter 


Values 


DETECT_MINAREA 


5 


DETECT.THRESH 


1.5(T 


ANALYSIS.THRESH 


1.5a 


FILTER_NAME 


tophat_3.0_3x3.conv 


DEBLENDJJTHRESH 


64 


DEBLEND_MINCONT 


0.001 


BACK.TYPE 


GLOBAL 


BACK_SIZE 


64 


BACK_FILTERSIZE 


3 


PHOT_APERTURES 


25 


BASIS.TYPE 


PIXEL_AUTO 


BASIS-NUMBER 


20 


PSF.SIZE 


25,25 



Table 2. Main input parameters set in the 
SExtractor and PSFEx configuration files. 
DETECT_MINAREA, BACK_SIZE, PHOT_APER.TUR.ES 
and PSF_SIZE are expressed in pixels. 
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18 2(1 22 24 2G 

Kron magnitude (B) 

Figure 1. /i max as a function of the Kron mag- 
nitude for stars (diagonal crosses) and galaxies 
(points) in the SExtractor catalog. 



5. Results 

In this section, we compare the results obtained 
using the two software packages. We will focus on 
four aspects of the extracted catalog namely: pho- 
tometric depth, reliability, accuracy of the derived 
photometry, and determination of the positions of 
the centroids. 

All the quantities and the statistics shown in 
this section are obtained by excluding saturated 
sources. In Fig. [TJ it is shown ^t max as a function 
of the Kron magnitude of the objects extracted by 
SExtractor. As evidenced by the flattening of star 
sequence, sources with magnitude B<19 mag are 
saturated in the simulated images. Starting from 
Tab.[3]and Fig. [2] we report the comparison among 
results obtained for the whole input magnitude 
range of unsaturated sources: 19-26 mag. How- 
ever, since we consider only input stellar sources 
recovered by both software packages and since the 
completeness limit of the DAOPHOT star catalog 
is B=24 mag (see Sect. 15. ll) . the last two reported 
magnitude bins are underpopulated and the re- 
sults may be affected by catalog incompleteness. 

5.1. Photometric depth 

The photometric limiting magnitude of the ex- 
tracted catalog is defined as the magnitude limit 
below which the completeness drops down to 90%, 
where the completeness is the ratio of number of 
detected sources, Nd otoctc( j, and number of input 
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Figure 2. Ratio between detected and input 
sources for different magnitude the bins. The dot- 
ted and the solid lines refer to SExtractor used 
with a Gaussian and a top-hat filter respectively, 
while the dashed line refers to values obtained with 
DAOPHOT. 

sources, N input . 

With DAOPHOT, the photometric depth depends 
mainly on the threshold applied, while for SEx- 
tractor, it depends also on the deblending of the 
sources, and on the filter used for the detection 
(see Sect. __]). As discussed in Sect 14.11 and I4.2[ in 
order to fix the thresholding and deblending pa- 
rameters, we performed several tests, by visually 
inspecting the extracted sources, and finally we 
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Figure 3. The left panel shows the ratio between the numbers of detected and input stars as a function of 
magnitude bins, as obtained by SExtractor (solid line) and by DAOPHOT (dashed line); in the right panel 
are plotted the same quantities, but for galaxies. 



fixed the values reported in Tab.[T]and[2j Then, we 
compared the results of source extraction obtained 
using two different filters: a Gaussian (dotted line 
in Fig. [5]) and a "top-hat" function (continuous 
line in Fig. [5]). As shown in Fig. [5J using SExtrac- 
tor with a top-hat filter, we can improve the detec- 
tion of faint sources. In this case, the depth of the 
catalog is ~ 25.0 mag. Hence, we refer to this filter 
in all the tests performed with SExtractor and re- 
ported below. Figure [2] shows also the percentage 
of extracted sources per magnitude bin obtained 
using DAOPHOT (dashed line). With this soft- 
ware package, the completeness drops rapidly to 
very low values for magnitudes fainter than B = 
22.0 mag. However, this comparison is mislead- 
ing. In fact, DAOPHOT is not designed to work 
with extended sources. For this reason, in Fig. [3a| 
we report the ratio between the detected sources, 
which are a priori known to be stars (S<j e tected), 
and the input stars (Si npu t)- In Fig. [3b] we show 
the same quantities but for galaxies (Gdetectedj 
Ginput)- We can see that the fraction of detected 
source is higher for stars for both SExtractor (B 
= 26.0 mag) and DAOPHOT (B = 24.0 mag). 
Hence, in conclusion, considering only stars, the 
final depth returned by DAOPHOT is - 2 mag 
brighter than those produced by SExtractor. 



5.2. Reliability of the catalog 

The reliability of the catalog is defined as the 
ratio between the number of well-classified sources 
and the number of the source s detected by the soft- 
ware packages (see Eq. (6) of lLaher. Grant fc Fuag 
2008). For these tests we use only the set of stars 
detected (Selected) and well-classified (Sdassffied) 
by both SExtractor and DAOPHOT. We com- 
pared results obtained with several methods to 
classify the sources. In fact, each method leads to 
a different estimate of the reliability. 
As far as DAOPHOT is concerned, we used the 
output parameters SHARP (see Fig. l4a]) and x ( see 
Fig. EE}, made the determination with ALLSTAR 
(see Sect. 13. lj) . Fig. l4al shows the distribution of 
ALLSTAR sharpness SHARP for our data. The 
separation between the two classes seems to be 
well defined. On the other hand, Fig. [4b] shows 
that the use of the \ parameter does not improve 
the star/galaxy classification. For this reason, we 
classified as stars all the sources with SHARP lower 
than 0. 

In order to investigate the reliability of the cat- 
alog made with SExtractor, we used both tra- 
ditional methods as well as the new parameter 
SPREAD_M0DEL. In Fig. EH we plot CLASS_STAR as 
a function of the Kron magnitude for our data. As 
shown, the lower the established limit to separate 
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Figure 4. Distribution of DAOPHOT sharpness (left panel) and x (right panel) as a function of the PSF 
magnitude for simulated stars (diagonal crosses) and galaxies (points). In the left panel the dashed line at 
zero sharpness is the adopted separation limit for the star/galaxy classification (see Sect. I5.2j) . 



stars and galaxies, the higher will be the contam- 
ination of the star subsample from galaxies. A 
reasonable limit for the separation is 0.98. 



SExtractor, the classifier SPREAD_M0DEL allows us 
to obtain a reliable star/galaxy separation down 
to B = 26 mag. 



Figures [5b] and [5c] show the locus of stars, 
selected according to the relation between half- 
light radius and /J, m ax, respectively, as a function 
of the Kron magnitude. There is an improve- 
ment of the source classification compared to the 
use of CLASSJ3TAR parameter, allowing a reliable 
star/galaxy separation down to B — 23.5 mag. 
Finally, Fig. [5d] shows SPREAD _M0DEL values as a 
function of Kron magnitude. Stars and galaxies 
tend to arrange themselves in two distinct places 
on the plot. Also in this case, the higher we choose 
the separation limit, the higher will be the con- 
tamination of the stellar sequence from galaxies. 
A good compromise between a reliable classifica- 
tion and a low contamination is the value 0.005. 
In Fig. [SJ it is shown the ratio between the sources 
correctly classified as stars using the stellarity in- 
dex (dotted line), spread model (continuous line) 
and sharpness parameter (dashed line), as func- 
tion of input magnitude. 

In conclusion, if we define a classification with a 
reliability of at least 90% with these methods we 
can acceptably classify the stars in DAOPHOT 
down to about 24 mag, which is the photometric 
depth of the extracted catalog, while in the case of 



5.3. Photometry 

In this section, we compare the results ob- 
tained with aperture and PSF photometry on 
the sample of stars detected by both SExtractor 
and DAOPHOT. We also investigate the results 
obtained with Kron, isophotal and model-fitting 
photometry for galaxies detected by SExtractor. 
In Tab. [3] we report the mean difference and the 
standard deviation between aperture and PSF 
magnitudes, as estimated by DAOPHOT (parts a 
and b, respectively) and SExtractor (parts c and 
d, respectively), against input magnitude. 
Figure [7] shows the residuals between aperture and 
input magnitudes (top panels), and the residuals 
between PSF and input magnitudes (bottom pan- 
els) , as estimated by DAOPHOT (left panels) and 
SExtractor (right panels). 

Tableland Fig. [7] show that there is a character- 
istic broadening of the residuals at fainter magni- 
tudes, as we expect when measurements become 
sky-noise dominated, but the spread in the case of 
PSF photometry remains smaller than for aper- 
ture measurements. This behavior is well known 



(e.g. iBecker et al.l 120071 ) for DAOPHOT, but it 
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Figure 5. Distribution of SExtractor stellarity index (panel a), half-light radius (panel b), /i ma x (panel c) 
and spread model (panel d), as a function of the Kron magnitudes for simulated stars (diagonal crosses) 
and galaxies (points). The dashed line in panels a and d is the adopted separation limit for the star/galaxy 
classification (see Sect. 15.21) . 
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is worth underlining that SExtractor has reached 
this level of accuracy in PSF photometry only af- 
ter the release of PSFEx. 

In the top part of the Tab. [3J we also report the 
mean difference and the standard deviation be- 
tween Kron (part a), isophotal (part b) and model 
magnitudes (part c), respectively, and input mag- 
nitudes for stars. 

For completeness, since SExtractor is designed 
also to obtain accurate galaxy photometry, we re- 
port in the bottom part of the Tab. [4] the mean dif- 
ference and the standard deviation between Kron 
(part d), isophotal (part e) and model magnitudes 
(part f), and input magnitudes for the "true" 
galaxies detected by the software package (See 
Sect. 5.2). 

By considering only stellar photometry, both soft- 
ware packages are able to deliver acceptable per- 
formances for both aperture and PSF photometry, 
up to a threshold two magnitudes brighter than 
the limiting magnitudes of the input simulated 
images, which is the completeness limit of the 
DAOPHOT catalog. Furthermore, the Kron mag- 
nitude yields ~ 94% o f the total source flux wit hin 
the adaptive aperture ( Bertin fc Arnoutslll996t ) so, 
accordingly, we see a shift of ~ 0.07 mag even in 
the brightest magnitude bin. On the other hand, 
the isophotal magnitude depends on the detection 
threshold and the model magnitudes (obtained 
through a sum of bulge plus disk), and produce 
an unbiased estimate of the total magnitude also 
for stars. 

In conclusion, the new PSF modeling of SExtrac- 
tor produces photometric measurements as accu- 
rate as those obtained with DAOPHOT. 



5.4. Centroids 

The last comparison is among extracted and 
input positions. There are different ways to ob- 
tain centroid measurements. As stated above, 
DAOPHOT can provide two different measure- 
ments for centroids. The simplest are the coor- 
dinates of the source barycenter, derived during 
the thresholding process. These coordinates can 
be redetermined by ALLSTAR, once DAOPHOT 
has built a PSF model, by applying a PSF correc- 
tion. 

Concerning SExtractor, we chose to compare the 
results obtained using the barycenter and the PSF 
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Input magnitude (B) 

Figure 6. Ratio between stars classified by Stel- 
larity Index (dotted line) and Spread Model (solid 
line) from SExtractor with threshold values re- 
spectively to 0.98 and 0.005, and by DAOPHOT 
sharpness (dashed line) with a threshold value 
equal to zero and detected stars, as function of 
input magnitude. 



corrected coordinates, as for DAOPHOT, and the 
results obtained by using the windowed positions 
along both axes. These coordinates are obtained 
by integrating pixel values within a circular Gaus- 
sian window. In Tab.0 it is reported the mean dif- 
ference between barycenter coordinates and PSF- 
corrected coordinates, estimated respectively with 
DAOPHOT (parts a and b) and SExtractor (parts 
c and d) and input coordinates. 
Finally, Tab. [6] shows the difference among input 
and windowed coordinates estimated by SExtrac- 
tor. 

Fig. [8] and [9] show the difference between the in- 
put and barycenter coordinates and between the 
input and PSF corrected coordinates. 
Both software packages show a bias between out- 
put centroid coordinates < 0.01 arcsec (equal to 
~ 0.47 pixel) and input X and Y, with an average 
deviation of < 0.02 arcsec (equal to ~ 0.94 pixel), 
down to the DAOPHOT completeness magnitude 
limit. These values are in particular improved in 
terms of average deviation (ctax(y) < 0.01 arc- 
sec), when PSF correction is applied. Hence, we 
can conclude that the results for centroids are sat- 
isfactory in both cases. 
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(c) (d) 

Figure 7. Top panels: Residuals between aperture magnitudes estimated by DAOPHOT (left panel) and by 
SExtractor (right panel), and input magnitudes for detected stars. Bottom panels: Residuals between PSF 
magnitude estimated by DAOPHOT (left panel) and by SExtractor (right panel), and input magnitude for 
detected stars. Superimposed red points and solid red lines draw the mean and standard deviation values 
reported in Tab. [3] 
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(a) (b) 

Figure 8. Difference between barycenter coordinates estimated by DAOPHOT (left panels) and by SExtractor 
(right panels), and input coordinates and as a function of input magnitude for detected stars. Superimposed 
red points and solid red lines draw the mean and standard deviation values reported in top part of Tab. [5] 
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(a) (b) 

Figure 9. Difference between PSF-corrected and input coordinates estimated by DAOPHOT (left panels) 
and by SExtractor (right panels), as a function of input magnitude for detected stars. Superimposed red 
points and solid red lines draw the mean and standard deviation values reported in bottom part of Tab. [5] 
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Bin 




Ani mC an 


CAm 


A III mean 


CAm 


Am m can 


f"Am 


Am m can 


f"Am 


(ma 


g) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(a) 


(b) 


(c) 


(d) 


19 - 


20 


0.007 


0.005 


0.002 


0.006 


-0.001 


0.003 


0.003 


0.003 


20 - 


21 


0.006 


0.010 


-0.001 


0.011 


-0.006 


0.005 


0.003 


0.004 


21 - 


22 


0.024 


0.070 


-0.011 


0.031 


-0.016 


0.011 


0.000 


0.012 


22 - 


23 


-0.020 


0.075 


-0.009 


0.059 


-0.023 


0.013 


0.003 


0.014 


23 - 


24 


-0.069 


0.157 


-0.011 


0.106 


-0.032 


0.026 


0.005 


0.025 


24 - 


25 


-0.147 


0.302 


-0.044 


0.250 


-0.034 


0.058 


-0.002 


0.055 


25 - 


26 


-0.306 


0.379 


-0.032 


0.444 


-0.122 


0.088 


-0.129 


0.129 



Table 3. The table reports, as a function of the magnitude bin (col. 1), the mean difference Am mcan 
(columns 2, 4, 6 and 8), and the standard deviation dAm (columns 3, 5, 7 and 9) between aperture and input 
magnitudes as estimated by DAOPHOT (part a) and by SExtractor (part b) . Parts c and d report the mean 
difference and the standard deviation between PSF and input magnitudes as obtained by using DAOPHOT 
and SExtractor, respectively. 



Bin 


Am mC an 


f"Am 


AlTlnican 


CAm 


A HI m ean 


C Am 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(a) 


(b) 


(c) 


19 - 20 


0.074 


0.005 


0.058 


0.007 


0.057 


0.023 


20 - 21 


0.077 


0.009 


0.073 


0.012 


0.048 


0.025 


21 - 22 


0.076 


0.027 


0.093 


0.024 


0.028 


0.046 


22 - 23 


0.076 


0.046 


0.128 


0.039 


-0.002 


0.074 


23 - 24 


0.087 


0.065 


0.216 


0.052 


-0.034 


0.091 


24 - 25 


0.051 


0.143 


0.389 


0.121 


-0.098 


0.146 


25 - 26 


0.147 


0.195 


0.749 


0.143 


-0.145 


0.151 



Table 4. The table reports, as a function of the magnitude bin (col. 1), the mean difference Am mcan 
(columns 2, 4 and 6), and the standard deviation a Am (columns 3, 5 and 7) between Kron (part a), isophotal 
(part b), model (part c) and and input magnitudes as obtained by using SExtractor. 
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Bin 




AX m can 


f"AX 


AY 

- 1 mean 


CTAY 


AX mean 


f"AX 


AY 

J- mean 


CAY 


(ma 


g) 


(pixel) 


(pixel) 


(pixel) 


(pixel) 


(pixel) 


(pixel) 




(pixel) 


(pixel) 


(a) 


(b) 


19 - 


20 


0.033 


0.055 


0.063 


0.062 


0.046 


0.011 




0.053 


0.016 


20 - 


21 


0.032 


0.045 


0.041 


0.059 


0.041 


0.019 




0.070 


0.054 


21 - 


22 


0.058 


0.061 


0.043 


0.050 


0.060 


0.024 




0.061 


0.024 


22 - 


23 


0.033 


0.068 


0.048 


0.058 


0.023 


0.072 




0.061 


0.047 


23 - 


24 


0.046 


0.108 


0.043 


0.090 


0.043 


0.062 




0.059 


0.079 


24 - 


25 


0.033 


0.209 


0.026 


0.202 


0.026 


0.132 




0.055 


0.138 


25 - 


26 


0.027 


0.167 


0.107 


0.288 


-0.044 


0.142 




0.124 


0.187 


(c) 


(d) 


19 - 


20 


0.050 


0.005 


0.050 


0.004 


0.051 


0.005 




0.050 


0.003 


20 - 


21 


0.052 


0.009 


0.051 


0.006 


0.051 


0.007 




0.051 


0.006 


21 - 


22 


0.052 


0.009 


0.053 


0.014 


0.054 


0.009 




0.053 


0.014 


22 - 


23 


0.044 


0.020 


0.054 


0.022 


0.044 


0.022 




0.052 


0.022 


23 - 


24 


0.043 


0.048 


0.046 


0.053 


0.044 


0.044 




0.049 


0.052 


24 - 


25 


0.031 


0.116 


0.043 


0.111 


0.033 


0.114 




0.047 


0.109 


25 - 


26 


-0.096 


0.149 


0.097 


0.189 


-0.007 


0.186 




0.107 


0.269 



Table 5. The table reports, as a function of the magnitude bin (col. 1), the mean difference between 
DAOPHOT X (col. 2),Y (col. 4) barycenter measure and input X,Y with the relative standard deviation 
(columns 3 and 5) in the part a, while in the part b there are the mean difference between SExtractor X 
(col. 6),Y (col. 8) barycenter measure and input X, Y with the relative standard deviation (columns 7 and 
9). In parts c and d arc reported the mean difference between X (columns 2 and 6), Y (columns 4 and 8) 
PSF corrected and input measurements obtained by using DAOPHOT and SExtractor, respectively, with 
the relative standard deviation (columns 3, 5, 7 and 9). 



Bin 




AXmcan 


0"AX 


AY 

1 mean 


CAY 


(ma 


g) 


(pixel) 


(pixel) 


(pixel) 


(pixel) 


19 - 


20 


0.050 


0.004 


0.050 


0.004 


20 - 


21 


0.051 


0.006 


0.051 


0.007 


21 - 


22 


0.054 


0.010 


0.056 


0.016 


22 - 


23 


0.032 


0.036 


0.054 


0.025 


23 - 


24 


0.044 


0.046 


0.052 


0.057 


24 - 


25 


0.028 


0.127 


0.040 


0.120 


25 - 


26 


-0.043 


0.124 


0.134 


0.173 



Table 6. The table reports, as a function of the magnitude bin (col. 1), the mean difference between X (col. 2) 
and Y (col. 4) windowed and input measurements as estimated by SExtractor with the relative standard 
deviation (columns 3 and 5). 
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6. Non uniform star distribution 

In order to evaluate the performances of both 
software in crowded fields, we tested them on a 
simulated image with a non uniform stellar distri- 
bution, i.e., showing an overdensity of stars in the 
center. In the left panel of Fig. QUI we plot the 
spatial distribution of the simulated stars, while 
the stellar density of the field, as a function of 
the distance from the center, is shown in the right 
panel. 

In order to carry out the test, we fixed the val- 
ues for DAOPHOT, ALLSTAR, SExtractor and 
PSFEx as in Sect.sOandO 
The obtained results were evaluated by following 
the same criteria as in Sect. [5] For all tests, we ex- 
cluded saturated stars, which are stars with mag- 
nitude brighter than B = 17 mag. 
Also, in case of a non uniform stellar field, we ob- 
tain that, for point-like sources, the final depth 
returned by DAOPHOT is ~ 2 mag brighter than 
those produced by SExtractor. In fact, SExtrac- 
tor can recover the input star catalog, resulting 
in a photometric depth of B = 25 mag, while 
for DAOPHOT the resulted limit in magnitude is 
around 23 mag. 

Moreover, by considering the reliability of the ex- 
tracted catalogs, it is possible to separate stars 
and galaxies down to the magnitude limit for each 
software package: B = 25 mag for SExtractor and 
B = 23 mag for DAOPHOT. 
In order to not bias the comparison of photometric 
measurements, we consider only the input stellar 
sources recovered by both software packages. As a 
first step, we analyze results for aperture photom- 
etry. Comparing the values of the mean and the 
standard deviation of the difference between aper- 
ture and input magnitudes with those reported in 
Tab. [H we obtained a decline in the quality of 
the results, in particular, those produced by SEx- 
tractor. In fact, in the magnitude bin 22-23, we 
go from Am ± <7A m — -0.009 ± 0.059 mag in the 
case of non-crowded field, to -0.165± 0.287 mag for 
SExtractor and from Am ± a Am -0.020 = 0.075 ± 
mag to -0.487 ± 0.474 for DAOPHOT. 
By using PSF-fitting photometry, there is an im- 
provement in the determination of magnitude with 
respect to aperture photometry, but the results 
show a quality decline as compared to the case of 
less crowded stars. 



By considering the last bin, for DAOPHOT PSF- 
fitting photometry, the results change from Am ± 
(J Am = -0.023 ± 0.013 mag to -0.079 ± 0.080 mag. 
While, by using SExtractor PSF-fitting photome- 
try, we pass from Am±<7Am = -0.003 ± 0.014 mag 
to -0.063 ± 0.110 mag. 

Finally, we compare the results among those ob- 
tained with barycenter and PSF corrected cen- 
troids and those reported in Tab. [5] In the last 
bin of magnitude, DAOPHOT measurements of 
PSF centroids go from AX ± a ax', AY ± a ay = 
-0.044 ± 0.020; -0.054 ± 0.022 pixel to AX±cr A x; 
AY ± a ay = -0.045 ± 0.146; -0.040 ± 0.148 pixel. 
On the other hand, in the same bin, SExtractor 
measurements of PSF centroids go from AX ± 
& ax i AY ± a ay = -0.044 ± 0.022; -0.052 ± 0.022 
pixel to AX ± ctax; AY ± a A Y = -0.063 ± 0.133; 
-0.055 ± 0.150 pixel. 

On the basis of this comparison, we can conclude 
that, as shown in Sect [SJ the results obtained by 
both packages, in the case of an overdense re- 
gion, are completely equivalent both for PSF fit- 
ting photometry and PSF fitting determination of 
centroids, and with regard of the reliability of the 
extracted catalog. With SExtractor, however, we 
can optimize the detection of the sources to re- 
cover the whole input catalog. 

7. Summary and conclusions 

The advent of new photometric surveys and the 
need to deal with fainter sources have led to an in- 
crease in the demand for quality and accuracy of 
photometric measurements. Moreover, the anal- 
ysis of massive data sets needs a single software 
tool in order to minimize the number of required 
reprocessing steps. Therefore, a crucial aspect in 
this context is the choice of source detection soft- 
ware package. An important innovation was in- 
troduced since 2010 with the public release of PS- 
FEx, which can work in combination with SEx- 
tractor. This new software package has deeply 
improved the performances of SExtractor, filling 
the gap against other available procedures, such 
as DAOPHOT and PHOTO, concerning PSF pho- 
tometry. 

In the present work, for the first time, the perfor- 
mances of DAOPHOT and SExtractor are com- 
pared to check the quality and accuracy of extrac- 
tion procedures with respect to the requirements 
of modern wide surveys. The software packages 
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Figure 10. Left panel: Spatial distribution of input stars for the over-dense stellar field. Colors of points 
refer to the surface number density of stars. Density values are reported in the vertical color-bar. Right 
panel: Numerical stellar density as a function of the distance from the center. Dotted line represents the 
average stellar density (~ 40 stars/arcmin 2 ) obtained excluding the circular region within a radius of 1.5 
arcmin from the center. 



are tested on two kinds of B-band simulated im- 
ages having, respectively, a uniform spatial distri- 
bution of sources and an overdensity in the center. 
In the first case, by considering only the number 
of extracted sources, it appears that the limiting 
magnitude for the extracted catalog is extremely 
low, in particular for DAOPHOT is B = 22 mag. 
If we limit to consider only stellar sources, the 
photometric depth is improved down to 24 mag 
for DAOPHOT and 26 mag for SExtractor. This 
could be related to the fact that only SExtractor 
gives the possibility to choose the filter for catalog 
extraction, according to the image characteristics 
(see Sect. 14. 2[) . In fact, as shown in Sect 15.11 the 
filter choice affects the photometric depth of the 
extracted catalog. 

A relevant aspect of the catalog extraction is the 
capability to discriminate between extended and 
point-like sources. As we have seen, within the dif- 
ferent software packages, there are various meth- 
ods to perform the star/galaxy classification. In 
particular the sharpness parameter available in 
DAOPHOT, improved by using ALLSTAR, re- 
turns a reliable star/galaxy classification down to 
the photometric depth of the catalog (B = 24 
mag). All the traditional methods available in 
SExtractor, instead, limit the star/galaxy classifi- 
cation, at least one magnitude above the complete- 



ness magnitude of the catalog. The new parame- 
ter SPREAD_M0DEL, which is a discriminant between 
the best fitting local PSF and a more extended 
model, has largely improved the classification, al- 
lowing to separate extended and point-like sources 
down to the completeness limit of the catalog (B 
= 26 mag for stellar sources). 
Since DAOPHOT is mainly designed to perform 
stellar photometry, in order to not bias the com- 
parison of photometric measurements, we consid- 
ered only input stellar sources recovered by both 
software packages. Both software tools are able to 
deliver acceptable performances in both aperture 
(with a (JAm < 0.2 mag) and PSF photometry 
(with a (JAm < 0.03 mag), down to B = 24 mag, 
the completeness limit of the DAOPHOT catalog. 
Moreover, since SExtractor allows us to derive dif- 
ferent estimates of the total magnitudes of sources, 
which we can also compare among themselves: 
Kron, isophotal and model magnitudes. The 
isophotal magnitude is highly dependent on the 
detection threshold. In fact, in the [23-24] magni- 
tude bin, there is a higher shift of Am (0.216 mag) 
than in other magnitudes. The Kron magnitude 
yields ~ 94% of t he total source flux with in the 
adaptive aperture ( Bertin fc Arnoutsl 19961) . Ac- 
cordingly, we find a shift of ~ 0.07 mag even in 
the brightest magnitude bin. The model magni- 
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tude results a good estimate of the input magni- 
tude also for stars, with an error of 0.091 mag in 
the [23 - 24] magnitude bin. 

An accurate determination of the object's cen- 
troids is crucial for the relative astrometry and 
thus, also for matching sources in different bands. 
Both software packages show a bias between out- 
put centroids and input X and Y coordinates < 
0.01 arcsec, with an average deviation of < 0.02 
arcsec down to B = 24 mag. These values are im- 
proved in terms of average deviation (<jax(y) < 
0.01 arcsec), when PSF correction is applied. So, 
we can conclude that these results are satisfactory 
in both cases. 

We also tested both software packages on simu- 
lated images with a non uniform stellar distribu- 
tion, i.e., showing an overdensity of stars in the 
center. By analyzing the results obtained in this 
last case, we can confirm the conclusions described 
above. In particular, DAOPHOT provides a cata- 
log ~ 2 mag shallower than the one extracted by 
SExtractor. On the other hand, by analyzing the 
extracted catalogs in terms of the mean difference 
and the standard deviation among output and in- 
put magnitudes and centroids, we notice a decline 
in the quality of the results for both software pack- 
ages, with respect to the case of a uniform spa- 
tial distribution of stars. Finally, DAOPHOT and 
ALLSTAR provide very accurate and reliable PSF 
photometry, with a robust star-galaxy separation. 
However, it is not useful for galaxy characteriza- 
tion. On the other hand SExtractor, associated 
with PSFEx, turns competitive in terms of PSF 
photometry. It returns acceptable aperture pho- 
tometry and accurate PSF modeling also for faint 
sources. The windowed centroids are as good as 
PSF centroids. Moreover, SExtractor allows to go 
very deep in source detection through a properly 
choice of image filtering masks. The deblending 
model is very extensible, and the use of neural 
networks for object classification, plus the novel 
SPREAD_M0DEL parameter, push down to the lim- 
iting magnitude the capability of star/galaxy sep- 
aration. Considering that SExtractor returns ac- 
curate photometry also for galaxies, we can con- 
clude that the new version of SExtractor, used 
in combination with PSFEx, represents a very 
powerful software tool for source extraction, with 
performances comparable to DAOPHOT also for 
overdense stellar fields. In the next years, it will 



be important to extensively test SExtractor plus 
PSFEx on real crowded stellar fields in order to 
definitively assess the performances of this soft- 
ware tool. However, an important aspect for the 
use of PSFEx and SExtractor, that we cannot 
avoid to mention, is the processing time. Without 
considering problems, such as degradation in per- 
formances during periods of heavy disk access, on 
average, SExtractor requires 0.5 s per detection to 
perform PSF photometry and source modeling, by 
using one single CPU with 6GB of RAM. This sug- 
gests that currently, the only disadvantage of us- 
ing SExtractor and PSFEx on wide field images is 
the processing time. However, on the other hand, 
although DAOPHOT is more efficient in terms of 
processing time just for the calculation, it requires 
more time if the user visually inspects the modeled 
PSF stars. 
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